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In the first part of this paper, an approximate theory is developed for the leading order concentra- 
tion dependence of the sedimentation coefficient for rod-like colloids/polymers/macromolecules. To 
first order in volume fraction tp of rods, the sedimentation coefficient is written as 1 + atp. For large 
aspect ratio L/D (L is the rod length, D it's thickness) a is found to very like oc (i) /log (tj)- 
This theoretical prediction is compared to experimental results. In the second part, experiments on 
/d-virus are described, both in the isotropic and nematic phase. First order in concentration results 
for this very long and thin (semi-flexible) rod are in agreement with the above theoretical prediction. 
Sedimentation profiles for the nematic phase show two sedimentation fronts. This result indicates 
that the nematic phase becomes unstable with the respect to isotropic phase during sedimentation. 

I. INTRODUCTION 

There is extensive literature concerned with sedimentation behaviour of spherically shaped colloidal particles (for 
a review see ref. [1]). Essentially exact predictions can be made for the sedimentation velocity of spherical colloids to 
first order in concentration [2]. For non-spherical colloids a similar exact prediction is non-existent. The only attempt 
to calculate the first order concentration dependence of the sedimentation velocity for rod like colloids we are aware of 
is due to Peterson [3]. This theory is based on approximate, orientationally pre-averaged hydrodynamic interactions 
between the colloidal rods and a rather crude estimate of certain multiple integrals that represent the ensemble 
averaged velocity. As yet there are no accurate expressions for hydrodynamic interaction tensors for rods. In the first 
part of the present paper, in section II, we calculate these interaction tensors in a mean-field approximation. In section 
III we use this approximate expression for the hydrodynamic interaction functions to derive an explicit expression for 
the first order in concentration coefficient of the sedimentation velocity as a function of the aspect ratio of the rods. 
This expression is found to agree remarkably well with Peterson's result for aspect ratios less than about 30. For larger 
aspect ratios our result for the first order in concentration coefficient is much larger then Peterson's prediction. In 
the second part of this paper, section IV, sedimentation experiments on /d-virus are discussed. Experiments are done 
at low concentration to find the first order concentration dependence, which is compared to the theory mentioned 
above. In addition, sedimentation experiments at larger concentrations, including the nematic phase are performed. 



II. HYDRODYNAMIC INTERACTION BETWEEN LONG AND THIN RODS 
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In order to calculate sedimentation velocities, the connection between translational and angular velocities, and 
hydrodynamic forces and torques must be found. In the present section such a relation will be established for two rods 
on the Rodne-Prager level, that is, with the neglect of reflection contributions between the rods. Considering only two 
rods limits the discussion on the sedimentation velocity to first order in concentration. Reflection contributions to the 
two-rod hydrodynamic interaction functions and multi-body rod interactions are both probably small in comparison 
to the Rodne-Prager terms, due to the fact that the distance between segments of different rods is of the order of the 
length of the rods, at least in the isotropic state. A Rodne-Prager approximation could therefore work quite well for 
long and thin rods, although explicit results for reflection contributions should be obtained to confirm this intuition. 

For the low Reynolds numbers under consideration, the translational velocities Vj, j = 1,2, and the angular 
velocities are linearly related to the hydrodynamic forces F /l and torques that the fluid exerts on the rods, 
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The superscripts "T" and "R" refer to translation and rotation, respectively, while the superscript "h" on the forces 
and torques refer to their hydrodynamic origin. On the Brownian time scale, the 3 x 3-dimensional mobility matrices 
M are functions of the positions of the centers of the two rods and their orientation. 

As will turn out, in order to find the sedimentation velocity, we need expressions for My T , for which approximations 
are obtained in subsection [IB. As a first step, the fluid flow field generated by a translating rod must be calculated. 
This is the subject of subsection II A. Rotation of rods also plays a role in sedimentation, but as will turn out, to 
first order in concentration and with the neglect of hydrodynamic reflection contributions, these do not contribute to 
the sedimentation velocity. Explicit expressions pertaining to the hydrodynamics of rotating rods are derived in the 
same spirit as for translating rods in appendix A. Subsection II C contains some concluding remarks. 

For the hydrodynamic calculations the rods will be thought of as a rigid string of spherical beads with diameter D. 
The length of the rods is L, and there are n + 1 = L/ D beads per rod, with n an even integer. 



A. Flow field generated by a translating rod 

The flow field generated by a rod that consists of n + 1 beads is given by, 
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u(r) = f dS'T(r-r')-fj(r'), (2) 



where T is the Oseen tensor, 
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with rjQ the shear viscosity of the solvent and f = r/r the unit position vector. Furthermore, fj(r') is the force per 
unit area that a surface element at r' of bead j exerts on the fluid, and dVj is the spherical surface of bead j. For long 
and thin rods, the distances r of interest, relative to the positions of the beads, are those for which r 3> D, with D 
the diameter of the beads. Now write r' = rj + R' with rj the position coordinate of the j th bead, so that R' = D/2, 
and Taylor expand the Oseen tensor in eq.(2) with respect to R'. Keeping only the first term in this Taylor expansion 
leads to relative errors of the order R'/r ~ D/L. Up to that order we then find, 
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u(r) = - £ T(r-rj).F]\ (4) 
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the total force that the fluid exerts on bead j. With the neglect of end-effects this force is equal for each bead, 
Fj = F h /(n + 1) = -^F' 1 , with F^ the total force on the rod. Eq.(4) thus reduces to, 

n n/2 

u(r) = - J2 T(r-rj).F h . (6) 

j=-n/2 

The force F h is calculated in terms of the translational velocity of the rod self-consistently from eq.(6) using Faxen's 
theorem for translational motion for each spherical bead, where the velocity vj of bead j is expressed in terms of the 
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force Fj on bead j and the velocity uo(rj) at the center of the bead that would have existed without that bead being 

where Vj is the gradient operator with respect to Vj. The first term on the right hand-side is just Stokes friction of a 
single bead in an unbounded fluid, while the second term accounts for hydrodynamic interaction between the beads. 
The fluid flow field Uq in turn is equal to, 
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where f * is the force per unit area that a surface element of bead i exerts on the fluid in the absence of bead j. For 
very long rods, consisting of many beads, the difference between fj (the corresponding force for the intact rod) and f* 
may be neglected: there are only a few neighbouring beads for which the difference is significant, but there are many 
more beads further away from bead j for which the difference is insignificant. To within the same approximations 
involved to arrive at eq.(4), eq.(8) can then be written as, 



n/2 

uo(rj) = - ]T T(rj-r,)-Ff. 
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Substitution of this expression into Faxen's theorem (7), and using that r,- 
of the rod, leads to, 
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Yi = (j — i)Du, with u the orientation 
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where eq.(3) has been used, together with, 
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For pure translational motion, the velocity vj of each bead is equal to the velocity v of the rod, so that both sides 
of eq.(10) can be summed over j, yielding for the left hand-side v L/D. Neglecting end-effects and replacing sums by 
integrals (which is allowed for long and thin rods), it is found that, 
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Notice that the Stokes friction contribution (the first term on the right hand-side in eq.(10)) is logarithmically small 
in comparison to the friction contribution due to hydrodynamic interaction between the beads. In fact, the Stokes 
contribution is neglected in eq.(12). A matrix inversion, in order to express F h in terms of v, and subsequent 
substitution into eq.(6), after rewriting the sum over beads as an integral over the center line of the rod, yields, 
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with r p the position coordinate of the rod. This is the approximate expression for the fluid flow generated by a 
translating long and thin rod that will be used in the following subsection to obtain an expression for the mobility 
matrices My T , j — 1,2. 



B. Calculation of M J 



In order to calculate the velocity v 2 that rod 2 acquires in the flow field (13) generated by a translating rod 

1, one should in principle perform a reflection calculation up to very high order : the field generated by rod 1 is 
scattered by each bead of rod 2 and subsequently reflected hence and forth between the different beads within rod 

2. Such a calculation is hardly feasible analytically. Here, the field generated by rod 1 that is incident on rod 2 is 
approximated by a constant fluid flow field u equal to the average of the incident field over the center line of rod 2. This 
"hydrodynamic mean-field approximation" is accurate for distances of the order L or larger, for which separations the 
incident field indeed becomes equal to a constant. For smaller distances between the rods this procedure provides a 
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semi-quantitative approximation. Within this approximation, the velocity of rod 2 immediately follows from eq.(12), 
with v = V2 — u, F h = F§ , the total force of the fluid on rod 2, and u = U2, the orientation of rod 2, 
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The average incident flow field follows from eqs.(13) and (12), with v = vi, the velocity of rod 1 and F h = F^, the 
force on rod 1, 
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where, 
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with T21 = Y2 — T\ the distance between the centers of the two rods. Notice that for distances Y21 between the centers 
of the rods larger than L, the matrix A asymptotes to T(r 2 i). By definition the following "mean-field" expressions 
for the translational mobility matrices are thus obtained (after an interchange of the indices 1 and 2), 
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One might try to device approximate expressions for the matrix A. However, sedimentation velocities are obtained 
as ensemble averages, also with respect to orientations, giving rise to integrals with respect to ri2 and ui,2, which can 
be evaluated by numerical integration. 



C. Concluding Remarks 

The approximations involved in the above discussion are justified for very long rods, since C(l)-constants are 
neglected against terms of order ln{L/D}, both by neglecting end-effects and replacing sums over beads by integrals 
(for the evaluation of the sums in Faxen's theorem in eq.(10)). Such approximations are most important for the 
diagonal mobility matrix (notice that factors \n{L/D} do not appear in the off-diagonal matrice Mj^\ due to 

the resubstitution of velocities in terms of forces). Both end-effects and the mathematical approximations involved in 
the calculation of the diagonal mobility matrix Mjr in eq.(17) may be accurately accounted for by the replacement, 
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with v = v±_ and v = v\\ a constant, pertaining to translational motion perpendicular and parallel to the rods 
orientation, respectively. This correction is experimentally significant for somewhat shorter rods (L/D < 20, say), 
but vanishes relatively to the logarithmic term for very long rods. The actual values of v± and for cylindrical rods 
are equal to [4], 



u± = -0.84, 
v\\ = 0.21. 
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A more accurate expression for than in eq.(17) is, 
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In the sequel we will use this expression for the mobility matrix MjJ 1 ' instead of eq.(17). The approximations involved 
in the off-diagonal mobility matrices in eqs. (34,38) are primarily due to the mean-field treatment of the incident flow 
field. It is probably a formidable task to improve on these expressions. 

The use of the more accurate expression (22) for the diagonal translational mobility matrix also circumvents the 
practical problem of calculating volume fractions of colloidal rod material from given values for L, D and number 
density. For the bead model it is not so clear how the volume of a rod must be expressed in terms of L and D. The 
volume of the cylindrical rod is simply equal to jD 2 L. 



III. AN EXPRESSION FOR THE SEDIMENTATION VELOCITY OF RODS 

The sedimentation of colloidal material induces, through the presence of the walls of the container, backflow of 
solvent. The backflow velocity is inhomogeneous, and varies on the length scale of the container. On a local scale, 
however, the backflow may be considered homogeneous, and the sedimentation velocity can be calculated relative 
to the local backflow velocity. This relative sedimentation velocity is a constant throughout the container (except 
possibly in a small region of extent L near the walls of the container, where gradients of the backflow velocity are 
large), and depends only on the properties of the suspension. A formal evaluation of the sedimentation velocity 
directly from eq.(l), by ensemble averaging, leads to spurious divergences, which are the result of the neglect of the 
hydrodynamic effects of the walls of the container which lead to solvent backflow. Batchelor was the first to deal with 
these divergences correctly, and we will use his arguments here [2]. 

Ensemble averaging of V! in eq.(l) gives the sedimentation velocity v s , which is thus found to be equal to, 

v s = - < M T ! T • + p VM^ 2 T • F£ + M™ • Tf + pVM™ • t£ > , (23) 

where the brackets < • • • > denote ensemble averaging with respect to positions and orientations of the rods. The 
factors pV — N « N— 1 account for the presence of TV— 1 rods which all interact with rod 1 under consideration. The 
divergence problems mentioned above arising in the explicit evaluation of the ensemble averages will be dealt with 
later. 

In order to be able to calculate these ensemble averages, the forces and torques must be expressed in terms of the 
positions and orientations of the rods. On the Brownian time scale there is a balance between all the forces and 
torques on each of the rods, that is, the total force and torque are equal to zero. The total force in turn is equal 
to the sum of the force F^ that the fluid exerts on the rod, the interaction force Fj = — Vj$ (with $ the total 
interaction energy of the rods), the Brownian force F^ r — — fceTVj ln{P} (with k B Boltzmann's constant, T the 
temperature and P the probability density function for positions and orientations), and the external force F ext due 
to the gravitational field. Hence, 

F£ = Vfi + k B TV 3 \n{P) - F ext . (24) 

Similarly, the total torque is the sum of the hydrodynamic torque T^, the interaction torque —lZj$>, the Brownian 
torque —ksTTZj ln{P}, while the torque on each rod due to the homogeneous external force vanishes. Hence, 

T£ = -R.fi + k B TTZ 3 \n{P} , (25) 

where the rotation operator is defined as, 

£;(••■) = fij x V„, (••■), (26) 

with V Uj . the gradient operator with respect to Uj. Substitution of eqs. (24,25) into eq.(23) for the sedimentation 
velocity yields, 

v s = - < M^ T • [Vi* + k B TVi ln{P} - F ext ] 

pVMj? • [V 2 * + k B TV 2 ln{P} - F ext ] (27) 
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The next step in the explicit evaluation of these ensemble averages is to determine the stationary probability density 
function P = P(r\, r 2 , Ui, U2) for the positions and orientations of two rods. At this point it is convenient to introduce 
the pair-correlation function g, defined as, 
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P(ri,r 2 ,ui,u 2 ) = P(ri,Ui)P(r 2 ,U2)g(ri2,Ui,u 2 ) , 



(28) 



where P(rj,Uj) is the probability density function for the position and orientation of a single rod. For spherical 
particles in a homogeneous external gravitational field, the probability density function for the position coordinates is 
simply the equilibrium function, without an external field, provided that the particles are identical. The probability 
density function differs from the equilibrium function only in case the relative sedimentation velocity of two spheres is 
different, for example due to differing masses and/or sizes. For rods, things are somewhat more complicated. Even if 
two rods are identical their relative sedimentation velocity generally differs as a result of the fact that the translational 
friction constant of rods is orientation dependent (see eq.(17)). The probability density function is generally dependent 
on the external force due to the fact that rods with different orientations overtake each other during sedimentation. 
Suppose, however, that the sedimentation velocity is so slow, that during a relative displacement of two rods in the 
gravitational field of the order of the length L of the rods, each rod rotated many times due to their Brownian motion. 
The relative sedimentation velocity of two rods then averages out to zero. For such a case, the pair-correlation function 
is only weakly perturbed by the external field, so that we may use its equilibrium Boltzmann form, 

#(ri,r 2 ,ui,u 2 ) = exp{-/3$(ri,r 2 ,ui,u 2 )} , (29) 

Let us derive the inequality that should be satisfied for eq.(29) to be valid. According to eq.(17), the largest relative 
sedimentation velocity Av s of two rods is approximately equal to, 

|Av.|« ^£>|F-| . (30) 
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On the other hand, the time T rot required for a rotational revolution is equal to, 
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Trot ~ ik B T\n{L/D} ■ 1 1 

The condition under which eq.(29) for the pair-correlation function is a good approximation is therefore, 
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> 1. (32) 



Hence, the work provided by the external force to displace a colloidal rod over a distance equal to its length should 
be much smaller than a few times the thermal energy of the rods. Substitution of typical numbers shows that 
this inequality is satisfied under normal practical circumstances.^] Furthermore, when alignment of the rods during 
sedimentation in a homogeneous suspension is of no importance, the one-particle probability density functions in 
eq.(28) are both constants equal to, 

P(r j ,u j ) = ^ v . (33) 

We will restrict ourselves here to the most common situation where the inequality (32) is satisfied, and assume 
negligible alignment, so that the probability density function is well approximated by eqs. (28,29,33). In that case 
many of the terms in eq.(27) for the sedimentation velocity cancel : the interaction contributions cancel against the 
Brownian terms. The sedimentation velocity reduces simply to, 

v s = [< M£ T > + p V < M?2 T >] • F ext . (34) 

Since Mj 2 T (ri 2 , Ui, U2) ~ l/i"i2 for large distances, its ensemble average with respect to position coordinates diverges. 
Such a spurious divergence is also found for spherical particles, and is the result of the neglect of the hydrodynamic 
effect of the walls of the container. Batchelor [2] showed that a formally divergent quantity, which is unambiguously 
finite valued on physical grounds, can be subtracted from the ensemble average, rendering a perfectly well defined 
sedimentation velocity. This finite valued quantity is formally divergent for the same reason that the ensemble average 



1 For example, for rods with a length of 100 nm, the sedimentation velocity should be much less than 1 mm/min, in order 
that the inequality (32) is satisfied. 
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of Mj^ is divergent, and subtraction accounts for the local hydrodynamic effects of the wallj^j Batchelor's argument 
is as follows. First define the velocity u(r | ri, • ■ • , rpj, ui, ■ • ■ , un) as the velocity at a point r (either in the fluid or 
inside a colloidal rod), given the positions n, • • -nv and orientations iii, • • ■ , ujv of N rods. In the laboratory reference 
frame the net flux of material through a cross sectional area must be zero. This means that the ensemble average of 
u must be zero. Hence, 

= < u(r |n, • • • ,r N ,ui, ■ ■ • ,u N ) > (35) 
= J dri ■ ■ J dr N <j> dui ■ ■ ■ <j> du N u(r|n, • • • ,r N ,fli, • • ■ , u N )P(n, , • • • , r N , Uj, • • ■ ,u N ) , 

where -P(n, • • ■ , u^v) is the probability density function for {n, ■ ■ • , ujv}- Formally, this ensemble average diverges for 
the same reason that the sedimentation velocity diverges : the flow field u varies like 1/r for large distances due to 
its Oseen contribution. The formally divergent expression (35), that must be zero for physical reasons, is subtracted 
from eq.(34) for the sedimentation velocity to render this expression convergent. The field u can generally be written 
as the sum of a two terms : a term to which the field would be equal to in the absence of reflection contributions plus 
a term that accounts for the reflection contributions. To within our approach reflection contributions are neglected 
so that only the former term survives here. Hence, 

N 

u(r |n, • • • , rjsr, ui, • • • , Ujv) = u(r — Tj) , for r in the fluid , 

.7=1 

= v s , for r in a core , (36) 

where the field u(r —Tj), for r in the fluid, is the sum of the flow fields in eq.(13), (with v replaced by v s — u s ), and 
eq.(66) in appendix A, (both considered as a function of the relative distance to the j th rod under consideration), 
that is, 

u(r - n) = u T (r - n) + u R (r - Tj) , (37) 
where uy(r — Tj) is the fluid flow generated by a translating rod as given in eq.(13), 

r L/2 
l-L/2 



u T (r - Tj ) = j f L/2 dlj T(r - Tj - Ijiij) • F ext , (38) 

L J-L/2 



where the inverse of eq.(12) is used, together with F h = —F ext , and u^(r — r 3 ) is the field generated by a rotating rod, 
which is similarly given by eq.(66) in appendix A. Operating on both sides of eq.(35) with J dr <f duP(r, u), where 
P(t, u) is the constant specified in eq.(33), and subtraction of the resulting equation from eq.(34) for the sedimentation 
velocity yields, for identical rods and to first order in concentration (rename ri = T2, r = n, Ui = U2 and u = iii), 

v s = -^v s + [< M^ T > +p V < Mj? >] • F ext 

- J dTl2 j diil j dil2 [ u r( r i2) + Ufl(r 12 )] X/(n I r 2 ,u 2 ) , (39) 

where Xf is the characteristic function that restricts the integrations to points r which are in the fluid, not inside the 
core of rod 2, 

X/(r|r 2 ,U2) = 1 , for r in the fluid , 

= , for r in the core of rod 2 . (40) 

Without interactions the angular velocity of each rod is simply proportional to the hydrodynamic torque on the same 
rod (see eqs. (59,64) in appendix A), which hydrodynamic torques are zero. Since the integral in the above expression 



2 To within the approximations made here, valid for long and thin rods, we do not encounter a conditionally divergent 
contribution due to terms ~ l/?"i2i as for spheres. Such a divergence can be dealt with by noting that the ensemble average 
of the deviatoric part of the stress tensor vanishes. We will not go into the extension of Batchelor's arguments to rods to deal 
with this conditionally divergence problem. 
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for the sedimentation velocity is multiplied by the concentration p, it follows that the rotational field u# does not 
contribute to first order in the density. For the same reason, in each term that is multiplied by the concentration, ~F ext 
may be expressed with eq.(22) in terms of the sedimentation velocity v° without interactions, at infinite dilution, as, 



v? = 
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We thus find the following expression for the sedimentation velocity, valid to first order in volume fraction, 
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where the functions /i and ji are equal to, 
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(44) 



where the expressions (13) and (38) for the mobility matrix M.J2 and the field u-r are used, respectively. We also 
used that integrals over the Oseen tensor must be proportional to the identity tensor, so that in these integrals the 
Oseen tensor may be replaced by the trace |Tr{T} of that tensor, which, according to its defining equation (3), is 
equal to, 



7>{T(r)} = 
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(45) 



For rods interacting only via a hard-core repulsion, it is shown in appendix B how to reduce the number of integrations, 
leading to the following results, 
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h = \^ + o(v/c) . 



(46) 
(47) 



The numerical value for /1 has been obtained by numerical integration and applies for hard-core interactions, where 
g is equal to when two cores overlap, and equal to 1 otherwise. The result for / 2 is independent of the kind of direct 
interaction between the rods. 

Substitution of the numerical values for /1 and f 2 from eqs. (46,47) and into eq.(42) gives our final result for the 
sedimentation velocity up to OiV/ C) contributions, 
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The volume fraction prefactor, 
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is plotted as a function of L/D in Fig.l (where both v±_ and v\\ are taken equal to 0). Also plotted is an older result 
due to Peterson [3] who predicted, 



8(3/8) 2 / 3 (L/J)V3 L 
a= 2HL/D} D- (50) 

In this latter theory back flow is not correctly accounted for, hydrodynamic interactions are orientationally preav- 
eraged and certain integrals are not precisely calculated but only estimated. The data points shown in Fig. la are 
experimental results for silica rods, coated with stearyl alcohol and disolved in cyclohexane. The data point x is 
taken from Ref. [5]. The point o is an unpublished result from the same author's. The data point if Fig. lb is a data 
point for fd- virus at high salt concentration, as obtained in the experimental section of the present paper. As can be 
seen from Fig. 1, the present prediction is virtually equal to that of Peterson for L/D < 30, but large differences are 
found for large aspect ratios. For large aspect ratios, a is predicted to vary like ~ (L/D) 2 /ln{L/D}, in contrast to 
Peterson's result ~ (L/D) 4 ^ 3 / ln{L/D}. For smaller aspect ratios a approaches approximately Batchelor's value for 
spheres a — 6.55, which is probably fortuitous in view of the approximations made here which limit the results to be 
meaningful only for long and thin rods. 



IV. EXPERIMENTAL RESULTS 

The concentration dependence of the sedimentation velocity predicted by Eq. (^) differs significantly from Peter- 
son's result Eq. ( [50"|) only for rods with large L/D. In our experiments we have used filamentous bacteriophage fd 
which is a rod-like virus with L/D ks 130. Other relevant physical characteristics of fd, are it's length L = 880nm, it's 
diameter D = 6.6 nm [6], and it's density of 1.285 mg/ml [7]. Because of it's large L/D ratio the virus is a semi-flexible 
rather then a rigid rod characterized with persistence length of 2.2 fim [8]. It's linear charge density is 10 e~/nm at 
pH8.2[9]. 

We have grown fd virus according to standard procedures of molecular biology described in Maniatis [10]. The virus 
suspension was first purified in a cesium chloride density gradient and then extensively dialyzed against tris buffer at 
pH 8.15 and at the desired ionic strength for the sedimentation experiments. After that the virus was concentrated by 
ultracentrifugation and from this stock solution a series of samples with different concentrations were prepared. The 
sedimentation velocity was measured on a Beckman XL-A analytical ultracentrifuge equipped with UV absorbance 
optics. Most of the experiments were done at 25 °C and at a centrifugal force equal to 45,500 g's (25,000 rpm). 
Before each sedimentation experiment the sample and rotor were allowed to equilibrate at the desired temperature 
for a few hours. Sedimentation data showed some unexpected features, interfering with straightforward calculation 
of the sedimentation coefficient. For this reason we have added appendix C, where the detailed analysis of our data 
is given. 

The measured sedimentation velocity for a range of volume fractions of fd from dilute solution up to a stable 
nematic phase are shown in Fig. ^|. All the samples in these measurements were kept at 8mM ionic strength. The 
sedimentation velocity of rods in the isotropic phase uniformly decreases with increasing concentration. After Ref. 
[10] we have tried to fit our experimental data to a functional form S$ = So(l — p<p) v where So is the sedimentation 
velocity at infinite dilution and <f> is the volume fraction of rods. The experimental values of sedimentation velocity 
are reported in Svedbergs where 15 = 10~ 13 s _1 . As seen from Fig. ||we obtain a reasonable fit to the experimental 
data in the isotropic phase and find that the sedimentation velocity at infinite dilution is Sq — 46.0 for the value of 
constants v = —1/3 and p=3600. After linearizing our fitted formula we find that the volume prefactor a « 1200 is 
much larger then predicted by Eq. (|49|). The reason for such a high value of slope a is the low ionic strength at which 
the experiment was performed. The same increase in a with decreasing ionic strength is observed in sedimentation of 
spherical particles [12,13]. Also we note that in this case the region where the sedimentation velocity varies linearly 
with rod concentration is limited to very low volume fraction of rods. 

It is a well known fact that elongated rods at high volume fractions undergo a first order phase transition to a 
liquid crystalline nematic phase [14]. The nematic phase is characterized by a short range liquid-like positional order 
and long range solid-like orientational order of rods, fd virus forms a cholesteric phase instead of the nematic phase 
[6,15]. Locally the cholesteric phase is equivalent to nematic, however on a macroscopic scale the average direction of 
molecules in a cholesteric phase forms a helix. The free energy difference between a cholesteric phase and a nematic 
phase is very small and although our experiments are performed on the cholestric phase only, we expect that our 
results are generic and would hold for a nematic phase of hard rods as well. 
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Bacteriophage fd in the cholesteric phase exhibits qualitatively new behavior when placed in centrifugal field. Instead 
of a single sedimenting boundary and single plateau we observe two boundaries with two plateau's sedimenting at 
different velocities as shown in Fig. |[ To confirm that this change in sedimentation behavior is indeed due to the 
formation of the nematic phase we have made a sample which is co-existing between the isotropic and nematic phase. 
After the sample had phase separated into macroscopically distinct co-existing phases, each phase was sedimcntcd 
separately. In the isotropic phase there was no sign of a second boundary, while in the nematic (cholesteric) phase we 
observed a fast sedimenting second boundary that was slightly more concentrated then the first component. On one 
hand, the slow component had a plateau concentration and a sedimentation velocity that was almost independent of 
the average concentration. On the other hand, the sedimentation velocity of the faster moving component rapidly 
decreases with increasing fd concentration and at the same time the difference between the plateau concentrations of 
the two components increases with increasing average concentration. 

The unstable sedimentation of colloidal rods in the nematic phase has a similar origin as the self-sharpening effect 
described in Appendix C. The reason for the instability is the discontinous jump in sedimentation velocity that 
occurs at the isotropic-nematic phase transition as is shown in Fig. ^ [16]. The denser nematic phase sediments at a 
significantly higher velocity then a more dilute isotropic phase. Initially a stable nematic phase occupies the whole 
sample length. When the centrifugal field is turned on a sharp sedimenting boundary starts moving towards the 
bottom of the container. Below this boundary (to the right) the rods are still in the nematic phase while above it 
the concentration of rods is very low and therefore they are in the isotropic state. This occurs because some particles 
will inevitably diffuse against the centrifugal field from a highly concentrated plateau into the dilute region. As this 
happens they simultaneously undergo a transition from the nematic to the isotropic phase. Since the sedimentation 
velocity of rods in the isotropic phase is much lower then in nematic phase, the probability of the rods diffusing from 
the isotropic phase back into the nematic phase is virtually zero. It is this asymmetry that results in a continuous 
flux of particles from the nematic into the isotropic phase and contributes to the formation of the second plateau in 
the isotropic phase that is moving at a slower speed. The concentration of the rods in the isotropic plateau will be 
very close to the concentration of isotropic rods co-existing with the nematic phase due to the self-sharpening effect 
because dilute isotropic rods will catch up with more concentrated isotropic rods. However, the highest concentration 
the isotropic rods can attain is the co-existence concentration between the isotropic and nematic phases, as long as 
the nematic phase sediments faster then the isotropic phase. Indeed, this is very close to what we observe in Fig. |^. 
Another experimental observation corroborating our explanation is that the sedimentation velocity and concentration 
of the slower isotropic plateau does not change significantly with average concentration of rods as seen in Fig |2|. 

Since the theory presented in this paper is valid only to first order in concentration of rods, to obtain an accurate 



value of the prefactor a in Eq. 49 we have made additional measurements in the dilute to semi-dilute range. Our 



results are presented in Fig. JlL We note that the overlap to semi-dilute concentration for fd with L — 0.88^m is 
at volume fraction of 5.9 • 10 . Unlike the previous measurements, we have done these measurements at high ionic 
strength where the behavior of charged rods is expected to approach the behavior of hard rods. Additionally at high 
ionic strength we expect the sedimentation velocity to have a linear dependence on volume fraction of rods up to 
higher values of volume fraction. The results for ionic strength of 50mM and lOOmM ionic strength are shown in 
figure [|. The volume prefactor in Eq. ( ^9|) at 50 mM ionic strength is a = 450 ± 40 and at lOOmM ionic strength 
a = 440 ± 60 . We have repeated the experiment at lOOmM ionic strength on a different analytical Beckman Xl-A 
ultracentrifuge and obtained the following result a — 490 ± 50. We conclude that a — 470 ± 50 which is the result 
plotted in Fig. lb. Since the values of the coefficient a do not change much with changing ionic strength from 50 
mM to lOOmM we conclude that the charged rods have approached the hard rod limit. Note that because of it's 
large L/D ratio fd is slightly flexible with a persistence length which is 2.5 times it's contour length [7]. Still for 
the experimentally determined parameters of fd, which are L = 880nm and D = 6.6nm, our experimental results 
compare favorably to the Eq. which predicts the value of a = 488 (see Fig. 1). In contrast, the previous result 
due to Peterson in Eq. B(J predicts a lower value of a = 288 
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APPENDIX A 



Flow field generated by a rotating rod 



Consider a rod with its center at the origin, which rotates with an angular velocity SI. The angular velocity is 
decomposed in its component perpendicular and parallel to the rods center line, 



Sin = uu • SI 



si 



(51) 
(52) 



Due to the linearity of the governing hydrodynamic equations, the flow fields generated by a rod rotating along Slj_ 
and may be calculated separately and added to obtain the flow field of the rod rotating along SI. 

Let us first consider a rod rotating with an angular velocity Slj_. The flow field that is generated by this rotating 
rod is given by the general equation (2). The relative change of the velocity of the beads is ~ For beads further 
away from the origin one may therefore consider the velocity over larger groups of beads as being virtually constant. 
The force on bead j is then proportional to its own velocity, 



F' 1 = -CSl± x rj = -CDjJ7_L x u. 



(53) 



where C is an as yet unknown proportionality constant. This expression is not valid for beads close to the center of 
the rod : for these beads the forces may have a different direction than their velocity. The fluid flow field generated 
by a long and thin rod, however, is primarily determined by the relatively large velocities of the beads further away 
from its center. We may therefore use eq.(53), except for relatively few beads close to the center and near the tips of 
the rod. Since r. 



jDii, the torque is thus found, to leading order in D/L, to be equal to 

n/2 



£ r iX Fj = -CD*± 



J=-n/2 



u x (Sl± x u) = -CD" 
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si 



(54) 



since Sl± is perpendicular to u. It is used here that Ylj=i J 2 ~ + 1) (2A + 1). First of all, the constant C is 

calculated self-consistently from Faxen's theorem in the form of eq.(10). Multiplying both sides of eq.(10) by rj-x, 
using that r 3 x Vj = j 2 D 2 Sl±, and summation over beads, leads to, 



-(-) dHi 

12 \D 



3tt?7oD 



CD 

87T?7o 



- g(L/D)0 : 



where the function g is defined as, 
9(L/D) = 



1 



n/2 n/2 

<n+ l) 3 

j = -n/2 i=~n/2,i^j 
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For long and thin rods the summations may be replaced by integrals, leading to, 

g{L/D) = \ \n{L/D} , 
6 



(55) 



(56) 



(57) 



up to leading order in D/L. Substitution of eq.(17) for the torque yields a single equation for C, yielding, again up 
to leading order, 



C 



Hence, from eq.(17), 



n 



ln{L/£>} 



31n{L/D} Th 



(58) 



(59) 



The flow field ui that is generated by a rotating rod may now be obtained from eq. (4) , to within the same approxi- 
mations that were discussed in the section [I A, as, 
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ra/2 . „ 2 n/2 

u ± (r) = - ^ T(r-r,)-F^ = £ T(r - r,) • (tl ± x j u) . (60) 

Replacing the sum over beads by a line integral, we thus find, 

47T7?n f L/2 

U± (r) = ln{L/D} J_ L/2 dl T(r ~ Vp ~ m) ' ( ^ X 1 fi) ' (61) 

where r p is the position coordinate of the rod. 

Next consider a rod rotating with an angular velocity Hy. For this case we have to resort to Faxen's theorem for 
rotational motion of a bead, which reads, 

^,| =-^3T J h + iv j xuo(r j ), (62) 

where, as in the translational Faxen's theorem (7), u is the fluid flow velocity that would have existed in the absence 
of bead j. The first term on the right hand-side is just Stokes rotational friction of a single bead in an unbounded 
fluid, while the second term accounts for hydrodynamic interaction between the beads. The important thing to note 
here is that the fluid flow generated by a single rotating bead is now equal to, 

u;(r) = (|^|) n||X(r-rj), (63) 

so that this fluid flow is along the entire center line of the rod. This implies that hydrodynamic interaction between 
the beads is unimportant for this case. For a long and thin rod rotating along its center line, each bead experiences a 
rotational friction that is practically equal to the Stokes friction, as if each bead where alone in an unbounded fluid. 
As a result, the total torque on the rod is simply the sum of the Stokesian torques on the beads, so that it follows 
immediately from Faxen's theorem (62) that, 

Furthermore, the total fluid flow uy is simply the sum of the fluid flows (63) generated by the rotating beads as if they 
were alone in an unbounded fluid, since hydrodynamic interaction between the beads is unimportant in the present 
case. Replacing the sum by a line integral thus yields, 

»,, W - X jC dl |r-r p -Zu|3 ("II X (r - • (65) 

The fluid flow u = ui + uy generated by a rotating rod with an arbitrary angular velocity ft = + 0|| follows 
by combining eqs.(51,52) and (60,65), 



47r?7o 



U « = RZM jC dl T(r - l *> ' ( " X 1 G) + X jC dl Ir-^-lup • x (r - r p)) . (66) 

This approximate expression will be used in the following paragraph to obtain an expression for the mobility matrices 
M™,j = l,2. 

Calculation of M TR 

In order to calculate the velocity v 2 that rod 2 acquires in the flow field (66) generated by a rotating rod 1, we 
apply, without further discussion, the same "mean-field" approach as in the previous section. The velocity v 2 is 
approximated by taking the fluid flow field generated by the rotating rod as a constant, equal to the average of the 
actual field over the center line of the rod. Hence, 



ln{L/D} 

V2 = U — ; 



I + U 2 U 2 



2 ) 



(67) 
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where the average flow field in terms of the torque on rod 1 follows from eqs.(66) and (59,64), with SI = Six, the 
angular velocity of rod 1 and T± the torque on rod 1, 

rL/2 r L/2 



u = 71 / dh f dh T(r 21 + l 2 u 2 - hux) ■ (h iix x T*] 

L J-L/2 J-L/2 



'-L/2 J-L/2 
I rL/2 ,L/2 1 

+ 5 72 / d/ i / d/ 2 1 — r~- T-^( r21 + h " 2 ) x • T i) ■ ( 68 ) 

8ttt] L 2 J_ L/2 J_ L/2 \r 21 + l 2 u 2 - hux \ 3 

By definition the following "mean-field" expression for the translational-rotational mobility matrices are thus obtained 
(after an interchange of the indices 1 and 2), ^ 

M™ = , (69) 

L/2 r L/2 



12 L4 



/ dh dl 2 l 2 u 2 x T(r i2 + hui - l 2 u 2 ) 

J-L/2 J-L/2 



L/2 J-L/2 
| [L/2 .L/2 1 

dh / dl 2 ■ ri — : j-^—^ [u 2 x (n 2 + h ui)] u 2 ■ (70) 

/2 J-L/2 



SnrjoL 2 J_ L/2 J_ L/2 \ r 12 + /i Ui - Z 2 u 2 | 3 



A non-zero contribution to Mj^ stems entirely from reflection contributions, since a pure rotation of a single rod 
in an unbounded fluid does not induce a translational velocity of the same rod. As mentioned before, reflection 
contributions are small in the isotropic state, since the typical distance between the beads of different rods is of the 
order L. 



APPENDIX B 

As a first step in the evaluation of the integrals in eq.(43) for fx, the Fourier transform of the Oseen tensor 
(T(k) = I — kk , with k = k/k) is substituted, and the integrations with respect to h and l 2 are performed, 

with the result, 

fx = ~ 87t 5 DL J dkk~ 2 <j) diix <j) du 2 J dr x2 [g(ri2, ui, u 2 ) - Xf ( r i I r 2, u 2 )] 

x cxp{ik • r 12 } jo(^Lk • Ui) j (^Lk • u 2 ) , (71) 

where, 

sinja;} 

jo(x) = . (72) 

x 

Consider the integral with respect to Y\ 2l 

1 = J drx2 [g(ri2,ui,u 2 ) - x/(ri |r 2 , u 2 )] exp{ik- r 12 } . (73) 

Replace the expression in the square brackets by (g — 1) + (1 — x/)- The integral over 1 — Xf is easily found to be 
equal to, 

/7T 1 
dr 12 [1 -x/(ri|r 2 ,u 2 )] exp{ik ■ r 12 } = -D 2 L j (-Lk • u 2 ) , (74) 

while the integral over g — 1 is equal to, 



3 The outer product a x A of a vector a and a matrix A is defined as the matrix with column vectors equal to the outer 
product of a and the column vectors of A. The outer product is thus taken with respect to the first index on A. 
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dri2 [ff(ri2,Ui,u 2 ) - 1] exp{ik-ri 2 } = -2DL 2 |u x x u 2 | j (-Lk • Ui)j (-Lk • u 2 ) 



(75) 



These results are valid for kD <C 1 (say kD < 0.2), while, in addition, eq.(75) is valid for orientations where 
^ <C| iii x u 2 |- As will turn out, the fcD-dependence is of no importance for long and thin rods, since convergence 
of the wavevector integral is assured by the fcL-dependent functions, which tend to zero for wavevectors for which, 
indeed, kD <C 1. Moreover, the angular integration range, pertaining to orientations where j-/ | iii x u 2 | is not 
small, vanishes for long and thin rods. Substitution of the results (74,75) into eq.(71) for fx, and noting that after 
integration over orientations the dependence on the direction k of the wavevector is lost, so that its direction may be 
chosen along the z-direction, yields (with x — hkL), 



ft 



f dx <j> d\\.i <j> du.2 j 2 (x z 2 ) 



|ui x u 2 | Jo(xzi) 



7T D 



(76) 



with Zj, j = 1, 2, is the z-component of Uj. The second term between the square brackets is an 0(D/£) contribution 
as compared to the first term and may be neglected. Transforming the orientational integrals to spherical coordinates, 
for which Zj = cos{Qj}, and using that (with \I/ = tpi — (^ 2 ), 



ui x u 2 |= 1 - (cos{9i} cos{8 2 } + sin{8i} sin{9 2 } cos{*})' 



1/2 



(77) 



finally yields eq.(46) for fx. 

Next consider the evaluation of the integrals in eq.(44) for / 2 . That the integrals are convergent follows from the 
Taylor expansion, 



ri2 - a 



1 

rx2 



a - V- 



ri 2 



1 1 

-aa : VV 

2 r 12 



(78) 



Using this in eq. (44) and integration with respect to Ui shows that the integrand varies like 
V 2 r^" 2 1 = for r\2 0. Following the same procedure as above one finds, 
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for large r 12 , since 



h = 



2ir s DL 



dk 



1 



x jo(^Lkz 2 ) 



i ,i 
dzi / dz 2 

-l J-i 

jo(hLkzx) - 1 



(27r) 3 5(k)- jD 2 Lj (iLkz 2 ) 



(79) 



where S is the delta distribution. The second term in the square brackets is easily seen to be 0(T>/C), using the 
same integration tricks as above for the evaluation of fx- For the evaluation of the delta distribution contribution, 
the integrand can be expanded in a power series expansion in fc. Using that jo(x) = 1 — x 2 /6 + • • •, results in eq.(47) 
for f 2 . 



APPENDIX C 



An analytical centrifuge measures the concentration of sedimenting colloid along the centrifugal field. From a 
single run in an analytical ultracentrifuge we obtain a time sequence of plots usually taken every few minutes. A 
representative sequence of these plots is shown in Fig. ^. Each plot in the series indicates the fd concentration as a 
function of radial position in the cell at that particular time. The concentrations of the dilute virus solutions were 
determined with the extinction coefficients of 3.84 mg _1 cm 2 at 270 nm [7]. For samples with higher concentration 
the solution is optically opaque at 270 nm and therefore we measure absorbance at progressively higher wavelengths, 
which corresponds to a lower extinction coefficient of fd. The sedimenting particles in Fig. ^| move from left to right. 
The water- air interface is indicated by a sharp peak located at radial position of 5.95 cm that is due to refraction 
by the meniscus. Note that this peak does not move as a function of time indicating that the container does not 
leak. As the rods start sedimenting towards the cell bottom, the region at the top of the solution (to the right of the 
air-water interface and to the left of the sedimentation front in Fig. |^) is depleted of virus as indicated by absence 
of absorbtion. Also the value of the concentration of rods in the plateau region, always to the right of the depleted 
region, is decreasing as the bulk of the sample moves towards the bottom of the container. The reason for this is 
that the walls of the cell are not parallel to each other, but instead follow the lines of centrifugal field in order to 
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minimize convective disturbances, an effect refered to as "radial dilution" [15]. Between the flat plateau region and 
the depleted region there is a sharp boundary. 

At higher concentrations of fd we observed the appearance of a sharp peak at the sedimenting boundary as is shown 
in Fig. |(]. The peak height increases with increasing concentration while the magnitude of the peak is independent 
of the wavelength and thus this peak cannot be due to absorption of the fd, which is wavelength dependent. The 
probable cause of the peak is the refraction of incident light due to the steep gradient in the virus concentration and 
hence the refractive index at the sedimenting boundary. As the incident light is refracted away from the detector, 
less light is collected by it and this results in apparent increased absorption of the sample. The peak at the water/air 
meniscus has the same origin. 

Two factors that determine the shape of the sedimenting boundary are the diffusion constant and the self-sharpening 
effect [17]. The diffusion of the particles leads to gradual spreading of the initially very sharp boundary. This 
diffusion of particles is countered by the self-sharpening effect, which is due to the concentration dependence of the 
sedimentation velocity. On one hand, any molecule lagging behind the boundary is in a more dilute environment 
and will therefore sediment at an enhanced velocity. On the other hand, the particles in the plateau region are in a 
more concentrated environment and their sedimentation will be retarded. As a consequence the boundaries will self- 
sharpen. In a suspension of elongated particles the self-sharpening effect will be much stronger then in a suspension 
of globular particles because the volume prefactor a in Eq. is much larger for elongated particles then for globular 
particles. The pronounced self-sharpening effect leads to hyper-sharp boundaries, which result in a steep gradient of 
refractive index which in turn causes the artifacts shown in Fig. ^. In globular colloids these effects arc usually not 
observed. 

In sedimentation analysis it is assumed that the rate of movement of the sedimentation boundary is approximately 
equivalent to the sedimentation velocity of the particles in the plateau (bulk) region. To compare results from different 
runs it is usual to express the sedimentation velocity in units independent of centrifugal force as follows: 

1 dr _ 1 rflnr 

S -—rd- t -^— (80) 

The sedimentation velocity unit is called a Svedberg (S), with 1 S = [10 3 sec~l]. We define r as the radial position 
at the sedimentation boundary where the virus concentration is equal to half the concentration of the plateau region. 
This quantity is easily obtained from experimental data for samples at low concentration. For samples at higher 
concentration, where we observe a peak at the sedimenting boundary due to refraction of light, we define r as 
the radial position of the highest point of the peak. A typical plot of the logarithm of r against uj 2 t used in the 
determination of the sedimentation constant is shown in Fig. ^. Surprisingly, we found out that a linear function 
provided an inadequate fit to our data. When the sedimentation data are collected between radial positions of 6.1 cm 
and 6.8 cm a polynomial of second order fits the data much better: 

In r = A + Bu! 2 t + CujH 2 (81) 

We introduce the experimentally observed sedimentation constant S r by combining Eq. ^l] and Eq. ^Tj| 

S r = ^f- = B + 2Cu?t (82) 
auj z t 

From this equation we see that the experimentally measured sedimentation velocity is not a constant but depends 
on the position of the measurement r of equivalently time uj 2 t at which the sedimentation front is found at position 
r. The reason for this unexpected behavior is not clear, but we assume it is an instrumental artifact. There is no 
physical reason to believe that the sedimentation velocity is a function of time or of radial position in the cell. In 
table 1 we see that the coefficient C, obtained when the quadratic polynomial in Eq. |8l] is fitted to data in Fig. ^a is 
independent of concentration. This is another indication that this artifact is due to the instrument. 

Theoretically, the constant B in eq. (j82[ ) should be equal to the concentration dependent Svedberg constant: 

S = #o(l-a<£) (83) 

and the constant C should be zero. Instead, we have found that the experimental Svedberg (eq. [8^) is described by: 

S r = 5 + offset = So + offset - S a^ (84) 

where "offset" depends on position in the centrifuge, but is independent of colloid concentration. So is the Svedberg 
constant of the rods in the limit of zero concentration. However, the value of slope 
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(85) 



is independent of radial position (or equivalently uj 2 t) where we evaluate Eq. |82| as is shown in table [I]. The measure- 
ment artifact only introduces a position dependent offset in the sedimentation velocity which affects the measured 
value of So- From a few measurements where we did not observe measurement artifacts (C = 0) we obtained the value 
of So = 47. Since this is in good agreement with previous measurements we use this value throughout our analysis 
[18]- 

It is important to note that the dependence of S r on position r shown in Fig. [7] is not due to the decreasing 
concentration of rods in plateau, which in turn is due to radial dilution. To show this we have made two measurements. 
In a first measurement we evaluated the sedimentation velocity at the point where the sedimenting boundary is close 
to the bottom of the container. At this time, due to radial dilution, the plateau concentration is about 70% of the 
initial concentration. In the second run our initial concentration was 70% of the concentration of rods in the first run. 
In this run we evaluated the sedimentation velocity right at the beginning of the run. We find that sedimentation 
velocities obtained in these two ways are vastly different, which indicates that the systematic errors described are not 
due to radial dilution. 
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FIG. 1. a) The dependence of the coefficient a on L/D ratio according to Eq. (49) and to Peterson [3]. The coefficient a 
is the first order correction to sedimentation velocity due to finite concentration of the colloidal rods. The two data points are 
for silica rods, b) The coefficient a for L/D < 200. The data point is for /d-virus, as obtained in the experimental section of 
the present paper 
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FIG. 2. Dependence of sedimentation velocity of rods measured in Svedberg's (Eq. p^ ) on the average volume fraction of fd 
at an ionic strength of 8 mM, pH 8.2. The equation of the fitted curve is S(cf>) = 46.0(1 + 3600<?f>)~ 1/ ' 3 . A sample made with a 
volume fraction between the two dashed vertical lines is unstable and will spontaneously phase separate into an isotropic phase 
at volume fraction of 0.0081 and a nematic phase at volume fraction 0.0093. When rods are sedimented in the nematic phase of 
initial volume fraction as indicated in the nematic region of the plot, we observe two sedimenting boundaries with two different 
velocities and two different concentrations. The sedimentation velocity of the faster moving component is indicated with filled 
circles, while the sedimentation velocity of slower moving boundary is indicated with filled squares. The actual volume fraction 
of the slower moving component is approximately constant at 0.0081 implying that that component is in the isotropic phase. 
The concentration of the faster component increases with the average concentration and is always concentrated enough to be 
a nematic phase. 
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FIG. 3. A concentration profile of sedimenting fd virus in a nematic(cholesteric) phase taken at two different times. Instead 
of a single moving boundary and a single plateau we observe two moving boundaries and two plateaus. Increased absorbance 
at the bottom of the container is due to the accumulation of the virus particles. Peak "b" marks the fast sedimenting nematic 
boundary while peak "a" marks the slow sedimenting isotropic boundary. The two curves are offset for clarity. The concentration 
of the initially uniform nematic sample was 13 mg/ml. The concentration of the co-existing isotropic and nematic phases at 8 
mM ionic strength is 10.5 mg/ml and 12 mg/ml, respectively. 
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FIG. 4. a) Concentration dependent sedimentation velocity S r (Eg. p4[ ) for fd at 50mM ionic strength. The data is fitted 
to a linear function S r = 45.3 — 209800. The solid line is given by Eq. |85f which yields aSo- The overlap volume fraction for 
fd with L=880 nm is 5.9 ■ 10 -5 . b) Concentration dependent sedimentation velocity for fd at lOOmM ionic strength. The data 
is fitted to a linear function S r — 45.9 — 204500. 
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FIG. 5. Data obtained from an analytical centrifuge. A time series of the fd concentration as a function of radial position in 
the centrifuge taken at 6.5 min intervals with "a" the first scan and "f" the last. The steep step in concentration represents the 
sedimentation front, which moves away from the centrifuge rotation axis with time. In this particular case the centrifuge was 
spun at 25000 rpm and the centrifugal (sedimenting) field points from left to right. The sharp peak "x" at the radial position 
of 5.95 cm is due to refraction by the air-water meniscus. Radial dilution accounts for the diminishing plateau concentration 
with increasing time. 
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FIG. 6. A series of plots of fd concentration as a function of radial position at time intervals of approximately 12.7 min 
with "a" the first scan and "d" the last. The difference between this series and those in Fig. ^|is that here the concentration of 
fd is higher and the rotation speed was 20,000 rpm. The peak seen at the sedimenting boundary is an artifact of the detection 
system and is due to the refraction of light at a sharp step in the refractive index at the sedimenting boundary. Radial dilution 
lowers the plateau concentration with time. A similar peak occurs at the air/water meniscus "x", which is stationary. 
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TABLE I. The values of constants obtained from 


Eq. H being 


fitted to data from Fig. [ 


k. The second column indicates 



the value of the sedimentation velocity 5* r evaluated at the start of the sedimentation experiment at to t — or equivalently 



which is independent of concentration, 
cm, 



If wc 



r=6.1cm. The third column indicates the value of parameter C in Eq. 
evaluate Eq. 



S2| for the sedimentation velocity at the end of the sedimentation experiment r=6.8 cm, we obtain the values of 
the sedimentation velocity shown in the fourth column. Note that the value of the slope aSo (Eq. ^ij) does not depend on the 
radial position. The value of aSo from the data evaluated at r=6.1cm is 20,500 and at r=6.8 cm is 21,000. We use the value 
So = 47 to obtain a. 
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FIG. 7. Position of the sedimentation boundary plotted against u> 2 t. The circles represent measurements which were taken 
approximately every 2.5 min. The lines represent the second order polynomial fit to the data (Eq. 81). The plot with a larger 
slope corresponds to the sedimentation of fd virus in lOOmM ionic strength at a volume fraction of 7.75T0" 5 at 25,000 rpm. 
The other plot is for a higher volume fraction of fd equal to 2.63- 10~ 4 . 
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